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Abstract 

Zero bias conductance per spin of nanotube double junction (NTDJ) is inves- 
tigated theoretically using the tight binding model, unrestricted Hartree-Fock 
approximation and non-equilibrium Green's functions. NTDJ consists of two 
metallic nanotubes joined by a piece of semiconducting nanotube, with the 
transition between the nanotubes made up of sets of 5 and 7 member car- 
bon rings. A quantum well forms in the central semiconducting NT region, 
bounded by Schottky barriers. Spin current occurs when Coulomb interac- 
tions raise the spin degeneracy of resonant levels in the quantum well. As 
long as an appropriate semiconducting NT length is chosen, spin direction 
can be controlled by gate voltage, i.e., NTDJ functions as a nano spin filter. 
The combination of nanotechnology based on carbon nanotubes (NTs) [[l],[| and spintron- 
ics [H is thought to be extremely promising for future technological innovations. Tsukagoshi 
et al. have demonstrated efficient spin injection from cobalt electrodes in NTs. Q However, 
diffusion of magnetic atoms from the electrode into the NT results in an increase in un- 
controllable spin flips due to the large spin-orbit interactions of ferromagnetic atoms. Since 
long spin relaxation times are needed, spin niters composed of atoms with small spin-orbit 
interactions are desirable. In this Letter, a spin filter composed entirely of carbon atoms, 
arranged as a semiconducting NT between two metallic NTs is proposed. Schottky barriers 
at the interfaces form a quantum well within the semiconducting NT, while the metallic 
NTs function as leads. |5|,[| When Coulomb interactions lift spin degeneracy in the quantum 



well, resonant current becomes spin polarized. To the best of the Authors' knowledge, this 
is the first study into such resonant spin current. |7],§] A possible method of constructing 
the NT junction is to allow open ends of (3m, 0) zigzag NT and (3m ± 1,0) NT to react; 
the NTs then become connected by adjacent regions of 5 and 7 member carbon rings, while 
still maintaining an sp 2 network. P, |IT||TT| , |T2T] This structure is referred to as the NT double 
junction (NTDJ). In order to examine spin current in such devices, zero bias conductance 
per spin is calculated for (9,0)- (8,0)- (9,0) NTDJ. Results are closely related to the fourfold 



degeneracy of the central semiconducting NT. [|13 
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FIGURES 

FIG. 1. Atomic structure of nanotube double junction (NTDJ), showing an (8,0) NT connecting 
two infinitely long (9,0) NTs. Length of the (8,0) NT region, L is 9a, i.e., L/(3a) = 3 where 3a is 
the length of the unit cell of (9,0) NT. Note the 5 and 7 member carbon rings in the (8,0)/(9,0) 
transition region. The NTDJ is divided into regions C, s and d, with regions s and d approximated 
as periodic. Gate surface and NTDJ's axis are taken to be parallel to each other, with a separation 
distance of five times the radius of (9,0) NT. The length of (9,0) NT regions included in calculation 
region C was 19.5a, although for illustrative purposes only 7.5a is shown. Geometries of the 
transition regions were taken to be the same, with the position where 5 and 7 member rings 
aligned assumed to be nearest to the gate. 

FIG. 2. Potential energy, i.e. diagonal terms of the converged Hamiltonian H c , as a function of 
position across the left half of region C for an NTDJ with L = 13. Values are averaged over each 
zigzag row. 

Using the unrestricted Hartree-Fock approximation, both the Hamiltonian H a and the 
density matrix p a are dependent on spin a =f, |, and can be expressed using the tight- 
binding approximation limited to the it orbitals as fl4[ 



= hi - Ui,^ + s id e UiM, k - (i/2)) , (i) 

k,cr 

where t^j = — t < when % and j are nearest neighbors, = otherwise. NTDJ is divided 
into the two asymptotic 'lead' regions s, d, and a scattering region C, as shown in Figure 
1, with Hamiltonian matrix elements assigned to one of these regions. The lead regions 
are approximated as having a periodic structure based on the charges and Hamiltonian 
elements in the NT unit cells at the edges of region C, as indicated by rectangles in Fig.l. 
The charges contribute to diagonal terms of H c , while Hamiltonian elements determine H s d 
and r Sj rf. For simplicity, only nearest neighbor elements are treated for off diagonal elements 
of H s d and r s d . s region, d region and gate electrode voltages are referred to as V s , Vd and 
V g , respectively. Since the gate electrode is a planar perfect conductor, and V g = in our 



model, Coulomb interaction can be expressed as U^j = U[\fi — fj\) — U(\fi — Rj\) , where 
Rj is the image of fj. fl5j The 'bare' Coulomb kernel without a gate is assumed to be fll4" | 



U(r) = 2M/s[l + 16(r/a) 2 . (2) 

Voltages can be expressed in another way, V s = (V s — V^)/2, Vd = —V s and V g = — (V s + 
Vd)/2, since the difference between, not the absolute value of, voltages is of importance. In 
this letter, relative voltage representation is used and negative V g causing electrostatic hole 
doping is considered. 

FIG. 3. Conductance per spin as a function of gate voltage V g for L =13. e± * and e±\ are 
shown by dashed and solid lines, respectively. 

The density matrix p is calculated from the retarded Green's function G, self energy 
Ap = t°(E + i5 — Hp)' 1 ^ and its imaginary component T° = i(A° — A^*), giving 



G a = (E - H a c -A a s -A 
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where /fd = l/(exp(/3(£' — Ep)) + 1) and ksT — 1/(3 — O.Olt. To calculate zero bias 
conductance, Ep = eV s = eVd = \e\V g . Using Eqs.(|lD and (H), H and p can be calculated self- 
consistently for a given V g . When self-consistency has been achieved, zero bias conductance 



for each spin dl u /dV s d is calculated from |T7 



dVcd h 



/oo 
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For small values of \V g \, motion of holes into (8,0) NT is suppressed by the band-gap, 
resulting in more positively charged (9,0) NT than (8,0) NT and inducing Schottky barriers. 
Fig. 2 shows potential energy as a function of position along the NTDJ for different spins 
and different values of V g . It can be seen that as V g becomes more negative, the (8,0) NT 
region becomes electrostatically doped, lowering the Schottky barrier. Neighboring peaks 
and dips seen in the transition region between (9,0) NT and (8,0) NT are caused by an 



excess of electrons in the 5 member carbon rings, and a deficiency in the 7 member carbon 



rings pOylofl . Note that while the gate electrode covers both C and lead regions, and is 
much larger than a nanometer, it still exerts control over the atomistic spatial profile of the 
barrier. 

Conductance per spin as a function of gate voltage V g is shown in Figs. 3, 4 and 5 for 
L =13, 7 and 6 respectively, where L is the length of (8,0) NT in units of 3a. The first 
resonant conductance peak appears at the gate threshold, V g th, with only a small tunneling 
current present when V g > V g th- For smaller values of L, enhanced tunneling effects increase 
hole density in the (8,0) NT region, making V g th more negative. 

FIG. 4. Conductance per spin as a function of gate voltage V g for L =7. e±j and e±^ are 
shown by dashed and solid lines, respectively. When V g = — 2.14t, the spin polarization vanishes 
while the conductance peak appears. This peak is not due to e± >a but comes from the level just 
under e_ denoted by e'. e' is shown by the solid line with crosses only when V g < — 2.14t. 

Fig. 3 also shows the four eigenvalues of H c near the valence band edge, denoted by e^, 
where i = ± is K corner point degeneracy. These can be represented as 



40) 



,-) - K)n_ itU + J (i)i) rij _ CT , (5) 

where J, K, n^ a and ef} are the Coulomb integral, exchange integral, occupation number 
and eigenvalue with no Coulomb interaction, respectively. Unlike finite systems or single 
electron tunneling, @,|l^] rii^ can take continuous values over the range to 1. From Fig.3, 



it can be seen that large splitting between the levels of the different spin states only occurs 
when ej jCr crosses the Fermi level Ep = \e\V g , i.e., when occupation numbers change. This 
indicates that splitting is a function of Jn^ a and Kni ;a , not of e^°\ Since the transmission 
rate T(E) = Trfr^GT^G" 7 *] exhibits resonant peaks when E ~ e^, as \V g \ increases two 
spin down peaks appear first, followed by two spin up peaks. The pairing of peaks with 
common spin is due to Hund's rule, which holds when K > e+* — ef?. 



For a given L, barrier thickness D decreases as hole number N = A — Y^io n i,v increases, 
as seen in Fig. 2, broadening the energy of the resonant state, W . At the same time, Ju j) 
decreases due to an increase in spatial width of the quantum well, L — 2D. The net effect 
is an increase in W — J(ij) resulting in a decrease in spin polarization. To examine spin 
polarization effects, total energy E tot is approximated by 

/ (0) w W \ ^ 

i,a ^ I I J i 

- Kn +j(T n- 7 a + J (+ _)n +iCT n_ a > (6) 

o" <t,<t' 

Equation (|6]) is derived from 

AE tot = An iia {W{n i7a -l/2) + e i>a ) , (7) 

where the energy region — W/2 < E — e i (J < +W(n ija — 1/2) is occupied so that n ija changes 
at the energy E = e i)(7 + W{n iiC — 1/2). 

For < iV < 1, n +) f and n + | initially decrease, because the degeneracy of is slightly 
lifted by presence of pentagon-heptagon pairs making < e+ . Therefore = 1 — a 
,n +i i = 1 — N + a and n_,f = n_ ( j, = 1. Thus, E tot = (J(+,+) — W)(N — a)a + C, with 
C independent of a. For a given N, we define a min as the value of a which minimizes 
E tot . Thus, spin current occurs when a min = 0. Although a min = N is another solution, 
infinitesimal Zeeman energy favors a m i n = 0. Thus the condition allowing spin down 
current at the first resonant peak (N = 1) is W|jv=i < J( +i+ )|jv=i. This condition is only 
satisfied when L > 4; since when L is too small, large tunneling effects mean that W — J 
is positive, even without electrostatic doping. As L increases, the value of N over which 
polarization occurs also increases. This can be understood by noticing that increases in 
W — J(i,j) per N become smaller as L becomes larger, because increases in hole density in 
(8,0) NT also decrease. 

When L > 4 and 1 < iV < 2,(n_ t f, n +t f, Ji_ t j, w +> |) = (1, 1 — a, 1 — N 2 + oc,0), where 
N 2 = N — 1 and < e_ _j_ on account of Hund's rule. Then 
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7=( e W - e V U) + ^+^(-,-)-%,-))/2 • (9) 

When 7 is assumed to be positive and nearly constant, level configuration depends on 
the sign of W — J(+~) — (2^/N 2 ). For a negative sign, holes are only injected into spin 
down levels, i.e., a min = 0. Otherwise holes are injected into both spin levels, i.e., a min = 
N 2 /2 — "f/(W — </(+,-))• When L = 5 or 6, W — </(+,-) — (2^y/N 2 ) increases, reversing 
sign as N 2 increases from to 1. The increase in W can be seen in the broadening of 
the peak corresponding to (— , j) in Fig. 5. As a result, spin up peaks overlap the second 
spin down peak, largely countering spin up polarization. For L > 6, on the other hand, 
W — J(+,-) — (2j/N 2 ) remains negative, and so the second spin down peak separates from 
the spin up peaks. 

FIG. 5. Conductance per spin as a function of gate voltage V g for L =6. e±^ and e±^ are 
shown by dashed and solid lines, respectively. 

In addition to the lower limits of L under which spin polarization can occur, upper 
limits also exist. As L increases for fixed N, both thickness of the barrier D and quantum 
well width L — 2D increase. The former reduces W while the latter causes a decrease in 
the Coulomb integral J(ij)- When is less than either W or k B T, spin polarization 
disappears. Because W decreases for larger L, k B T becomes a more important factor in 
determining the upper limit of L than it was for the lower limit. Figure 3 shows that at 
room temperature, the upper limit is in excess of 39a. 

In this letter, it has been shown that the nanotube double junction (NTDJ) functions 
as a nano spin filter without the use of magnetic atoms. An infinitesimal magnetic field 
is necessary only to determine which spin becomes unoccupied first. Conditions for spin 
current are: (1) Fermi level is near the band gap edges at operating gate voltage V g ; (2) 
length of the semiconducting NT region, L, is larger than a particular threshold length. In 
fact, there are two threshold lengths, L\ and L 2 . For L > L 2) the direction of spin can be 
controlled by V g . For L 2 > L > L 1; however, spin polarization is limited to a single direction. 
Furthermore there is no spin current when L\ > L. In the system studied, these thresholds 
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are L\ = 4 and L2 = 6. Although the Coulomb interaction U(r) of eq.(|2|) was chosen to be 



slightly smaller than estimated in Ref. |TJj to achieve convergence of the self-consistent loop, 
the Author believes that the parameters used are realistic enough and results valid enough 
to prove that NTDJ is a useful device that combines nanotechnology and spintronics. 
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